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Abstract 

An exact solution of the Einstein equations for a Bianchi -I universe 
in the presence of dust, stiff matter and cosmological constant, gener- 
alising the well-known Heckmann-Schucking solution is presented. 
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1 Introduction 



In spite of the development of the numerical methods for integration of the 
Einstein equations and an extensive application of the qualitative theory of 
differential equations to analysis of the behavior of the cosmological models, 
the construction of exact solutions still represents rather an attractive task. 
One can underline two important classes of the cosmological solutions known 
now. One of them is the class of exact solutions of Friedmann-Robertson- 
Walker isotropic cosmological models, which constitute a basis for compari- 
son of theoretical predictions with observations pp. Another important ex- 
act cosmological solution is the anisotropic Kasner solution for the empty 
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Bianchi-I universe j2], whose importance for theoretical physics is connected 
with its role in the description of the oscillatory aapproach to the cosmolog- 
ical singularity [3], which in turn appears very promising topic for the study 
in the string and M-theory context [3]. 

The Heckmann - Schucking |3] anisotropic solution for the Bianchi-I uni- 
verse in the presence of the dust-like matter is of special interest because it 
constitutes some kind of bridge between these two types of the cosmological 
solutions: in the vicinity of the cosmological singularity it behaves as a Kas- 
ner universe, while at the later stage of the cosmological evolution it behaves 
as an isotropic flat Friedmann Universe. The Heckmann- Schucking solution 
has the following form: for the Bianchi-I universe with the metric 

ds 2 = dt 2 - a 2 (t)dx 2 - b 2 (t)dy 2 - c 2 (t)dz 2 (1) 
filled with dust whose equation of state is 

P = (2) 
the functions a(t), b(t) and c(t) are given by the formulae 

a(t) = a t p i(t + t ) 2 / 3 - p \ 
b(t) = b Q t p *(t + t ) 2 / 3 - p \ 

c(t) = Co t P3 (t + to) 2/3 - p3 , (3) 

where the exponents pi,P2 and p^ are the well-known Kasner exponents PHE] 
satisfying relations 

P1+P2+P3 = 1, (4) 
pI+pI+pI = 1- (5) 

It is easy to see that the solution (jHJ) is close to the Kasner solution when 
t <C to- 111 the limit t ^> t all the functions a(t), b(t) and c(t) are proportional 
to t 2//3 , i.e. their behavior coincides with that of the flat Friedmann universe 
filled with the dust. The energy density of the matter is given by the formula 

- E ° (6) 



t(t + 1 ) 



where the constant E ^ and does not depend on the choice of the Kasner 
exponents. 
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In this paper we generalise the Heckmann - Schucking solution for the 
case of the Bianchi - I universe filled with the mixture of three perfect flu- 
ids: dust, stiff matter and a cosmological constant. The presence of the 
cosmological constant in combination with other types of matter is of special 
interest, because the recent discovery of the phenomenon of the cosmic ac- 
celeration jB] makes the presence of the cosmological constant or some other 
exotic type of matter which mimics some basic features of the cosmological 
constant enavoidable for any realistic cosmological model. Another reason 
which makes the study of the models with the cosmological constant interest- 
ing is the fact that inflationary theories of a very early Universe |7] contain 
an effective cosmological constant providing a period of a quasi-exponential 
expansion at the beginning of the cosmological evolution. 

The plan of the paper is the following one: in second section we integrate 
the Einsten equations for the Bianchi - I model filled with the mixture of 
three perfect fluids mentioned above; in the third section we consider different 
limit cases for our solution and rewrite it in the form maximally close to the 
canonical form of the Heckmann- Schucking solution (J3j). 

2 Integration of the Einstein equations for a 
Bianchi - I cosmology in the presence of 
matter 

We shall look for solutions of the Einstein equations for the Bianchi - I 
universe with the metric (|1J) filled with the matter with an isotropic matter 
with the energy-momentum tensor which has the form 



It is convenient to represent the fucntions a(t), b(t) and c(t) in the following 
form: 



where R(t) is the conformal factor while the functions a(t) and (3{t) charac- 
terise the anisotropy of the model. The components of the Ricci tensor have 



T v = diag(e, -p, -p, -p). 



(7) 



a(t) 
b(t) 
c(t) 



R(t)exp(a(t)+P(t)), 
R(t)exp(a(t) -P(t)), 
R(t)exp(-2a(t)), 



(8) 
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the following form: 



(9) 



Rl = -(^+2^ + 3^(a + P)+a + p), (10) 
R 2 2 = -(^ + 2^ + 3^a-$) + a-(3y (11) 

(i2) 

Due to the isotropy of the energy-momentum tensor (JJJ) one has 

-^i = -^2 = -^f- ( 13 ) 

From this equation it is easy to get equations describing for the asymmetry 
functions a(t) and (3{t): 

R 

a + 3-a = 0, (14) 
R 

'(3 + 3^ = 0. (15) 

The form of Eqs. ((H)) and ()15j) coincides with that for the Klein - Gordon 
equation for the massless sclalar field on a Friedmann background, whose 
energy-momentum tensor coincides with the energy-momentum tensor of the 
stiff matter [S] with the equation of state 

p = e. (16) 

From Eqs. (|T4^) . ()15|) it follows immediately that 

a = w (17) 

P=%, (18) 

where «o and (3q are some positive constants. 

The 00 component of the Einstein equations has now the form 

R 2 . 2 (3 2 4ttG 

R 2 = a + T + — £ - (19) 
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We shall consider a universe filled with a mixture of three perfect fluids: 
dust, the stiff matter obeying the equation of state (jT5|) and the cosmo logical 
constant, whose equation of state is p = —e. Choosing convenient normal- 
isation of constants one can represent Eq. ()19|) for the universe filled with 
this mixture in the following form: 

R 2 2 $ 2 A M S 

_ = a . + <L +A + _ + _ (20) 

where A is the cosmological term, while the constants M and S characterise 
the quantity of the dust and of the stiff matter in the universe respectively. 

After substitution into Eq. (|20p the expressions for a and (3 from Eqs. 
(fT7jl and (fTHj) we come to the following equation for the conformal factor R(t) 

R 2 A M S 

W = A + W + W' (21) 

where 

S = S + a 2 + ^. (22) 

It is possible to integrate Eq. (J2~T|) explicitly and the result is the following 
one: 

R 3 (t) = ^-(coshSv/At - 1) + \J^-sinh3VAt. (23) 

The initial conditions for the solution fl23(l are chosen in such a way to provide 
the fulfilling of the relation R(0) = 0. Substituting the solution (J2llj) into Eqs. 
(fTTj) and (fTS|) one has after taking integrals: 

a( ^3vk ln b^= ' (24) 

*)=3^ in (^l5fc)' (25) 

In the expressions fJ24jl . (j2*H|) we have omitted the integration constants, in- 
clusion of which implies only the multiplication of factors a(t),b(t) and c(£) 
by some factors, that is not important for the flat (Bianchi - I) model. 
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Now, substituting the formulae (|23|). and (j23|) into Eq. (jSJ) one comes 
to the following cosmological solution: 




The solution presented above has rather a cumbersome form. In the next 
section we consider different limit cases and rewrite Eq. ()26j) in a more 
convenient Heckmann-Schucking-like form. 



3 Limit cases and the Heckmann-Schucking- 
like form of the anisotropic cosmological so- 
lution 

Let us consider the limit case A = in Eq. (J26j) . i.e. the case when the 
cosmological constant in absent. One has the following expressions: 

a(t)= Qm)V(* + *o) § " Pi , 
6(t)= (|m) V(t + t )3- p2 , 
c(t)= (^M)V(t + t )3- ps . (27) 
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Here we have introduced the parameters 



Pl=n 



3M ' 

1 _ a + p 
3 3^ ' 



(28) 



P2 = ~ 



3 3^ 

p 3 =~ + — JL (29) 



3 3VS 



o 

The parameter to plays the same role as in the Heckmann-Schucking solution 
@, but now it depends not only on the anisotropy parameters a and (3q but 
also on the quantity of the stiff matter in the universe(see Eq. ([22))) . It is easy 
to see that the parameters pi,p2 and p%, which can be called "quasi-Kasner" 
parameters satisfy the relations 

Pi +P2+P3 = 1, (30) 

pl+pi+pl = i-q\ (si) 

where 

2 2 S 2S , . 

Q = a5 = 32 

3 + fit + 5 350 

The parameter g reflects the presence of the stiff matter in a universe and 
was introduced in [H] . Apparently, when the stiff matter is absent q = and 
the solution (J27j) exactly coincides with the Heckmann-Schucking solution © 
for a Bianchi - I universe filled with dust. The presence of the stiff matter 
simply transforms the relation (J5J) for the Kasner exponents into Eq. ([31)1 . 

Now, it convenient to rewrite the solution (|26|) using the Heckmann - 
Schucking time parameter t Q and the quasi-Kasner exponents. We shall use 
also the parameter 

H = y/X, (33) 

which is nothing but a Hubble parameter for the DeSitter universe with a 
cosmological constant A. Now the solution (|2lj|) looks like 

, A (My ( 3H t \ p ^ ( . 3H ty* ( . 3H t 3H t 3H ty~^ 
a(t)=[- + — J ) ( Sinh — + — COsh — J 
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Kt) = 1^1 ( 1 + = ) ( smh fsinh ^ + ^ cosh W 




2 / \ 2 2 2 

r(7) = ( — I ( 1 + ^ ) ( smh^f fsinh^ + ^ cosh ^ 



The solution for the Bianchi - I cosmology in the presence of the cosmolog- 
ical constant, dust and stiff matter, presented in the Heckmann-Schucking- 
like form (J3"4"j) is more convenient for the comparison with other solutions. 
We have already considered the case A = which reproduces the Heckmann- 
Schucking solution (jHJ) up to change of the relation (J3TJ) between the Kasner 
exponents, which occurs due to the presence of the stiff matter. For complete- 
ness we shall write down also the solution for the case when dust component 
is absent. Taking the corresponding limit in Eq. ()34|) we easily come to 



(3H t 
tanh ^° 

i / 3H t 
b(t) = sinh(3# *)^ (tanh 



2 

c{t) = sinh(3i7 £)^ (tanh ' • (35) 



pi- 



We shall write down also asymptotics of the solution (J34)) for large and 
small times. At large times t — *> only the cosmological constant is relevant 
and a complete isotropisation takes place. The corresponding asymptotics 
look like 

a(t) ~ exip(Hot), b(t) ~ exp(H t), c{t) ~ exp(H t). (36) 

At the beginning of the cosmological evolution when t — > only anisotropy 
and the presence of stiff matter determine the behavior of functions a(t), b(t) 
and c(t) . Their form is now 

a(t) ~ t p \ b{t) ~ i p \ c(t) ~ t P3 , (37) 

where the exponents pi,P2,Pz satisfy as usual the relations (|3*T|). (|3*T|). Let 
us notice that in the model under consideration there are two time scales: to 
and Hq 1 . One can write down also "intermediate" asymptotics. In the case 
when H Q t < 1 the solution for the time scale 

to < t < — (38) 
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looks like 

a(t) ~ i§, b{t) ~ t§, c(t) ~ t§ (39) 

and corresponds to the flat Friedmann universe filled with dust. For the case 
when H t > 1 and 

7T < * < *o (40) 

the solution coincides with that from Eq. (|36p. i.e. describes the flat DeSitter 
universe. 

In conclusion, notice that recently P| exact solutions for Bianchi - I uni- 
verses filled by stiff matter or in the presence of the cosmological constant 
were considered with use of the original technique |H| using analogies between 
Einstein and Newton cosmologies. The qualitative theory for the Bianchi -I 
universes filled with different types of matter was developed in [ID] . 

This work was supported by RFBR via grants No 02-02-16817 and 00-15- 
96699 and by International Scientific Program " Relativistic Astrophysics" by 
the Ministry of Industry, Science and Technology of the Russian Federation. 
The authors are grateful to V.A. Belinsky for useful discussion. 



References 

[1] L.D. Landau and E.M. Lifshitz , The Classical Theory of Fields (Perg- 
amon Press), 1979. 

[2] E. Kasner, Amer. J. Math. 43 (1921) 217. 

[3] V.A. Belinsky and I.M. Khalatnikov, Sov. Phys. JETP 29 (1969) 911; 
I.M. Khalatnikov and E.M. Lifshitz, Phys. Rev. Lett. 24 (1970) 76; V.A. 
Belinsky, I.M. Khalatnikov and E.M. Lifshitz, Adv. Phys. 19 (1970) 525; 
31 (1982) 639; C. W. Misner, Phys. Rev. Lett. 22 (1969) 1071. 

[4] T. Damour and M. Henneaux, Phys. Lett. B 488 (2000) 108; Phys. Rev. 
Lett. 86 (2001) 4749; T. Damour, M. Henneaux, B. Julia and H. Nicolai, 
Phys. Lett. B 509 (2001) 323. 

[5] O. Heckmann and E. Schucking, Newtonsche und Einsteinsche Kosmolo- 
gie, Handbuch der Physik, v. 53 (1959) 489-519. 

[6] S.J. Perlmutter et al, Astroph. J. 517 (1999) 565; A. Riess et al, Astron. 
J. 116 (1998) 1009. 



9 



[7] A. A. Starobinsky, Stochastic De Sitter (inflationary) Stage in the Early 
Universe, in Field Theory, Quantum Gravity and Strings, (Eds. H.J. De 
Vega and N. Sanchez, Springer- Verlag, Berlin) (1986) 107; A.D. Linde 
A D Particle Physics and Inflationary Cosmology (Harward Academic 
Publishers, New York), 1990. 

[8] V.A. Belinsky and I.M. Khalatnikov, Sov. Phys. JETP 36 (1973) 591. 

[9] P. Nurowski, E. Schucking and A. Trautman, Relativistic gravitational 
fields with close Newtonian analogs, in Einstein's Path, Essays in honor 
of Engelbert Schucking, ed. A. Harvey, Chapter 23, (Springer, New York, 
1999). 

[10] I.M. Khalatnikov and V.A. Belinsky, Qualitative Cosmology, in Prob- 
lems of Modern Experimental and Theoretical Physics, dedicated to 
Yu.B. Khariton, (Ed. A. P. Alexandrov, Nauka, Leningrad, 1984, in Rus- 
sian), 299-322; I.M. Khalatnikov, Qualitative Cosmology, in Current 
topics in Astro fundamental Physics, 1 course, Erice, Italy, 1-8 Sept. 
1991, (Eds. N. Sanchez and A. Zichichi, World Scientific, Singapore, 
1992), 47-66. 



10 



